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Abstract 



The electronic response of a composite consisting of aligned metallic cylin- 
ders in vacuum is investigated, on the basis of photonic band structure calcula- 
tions. The effective long-wavelength dielectric response function is computed, 
as a function of the filling fraction. A spectral representation of the effective 
response is considered, and the surface mode strengths and positions are ana- 
lyzed. The range of validity of a Maxwell-Garnett-like approach is discussed, 
and the impact of our results on absorption spectra and electron energy-loss 
phenomena is addressed. 
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I. INTRODUCTION 



Both theoretical and experimental investigations of the optical properties of composite 
materials have been of basic importance over the years!. In particular, effective medium 
theories, in which the electromagnetic response of the composite is described in terms of an 
effective dielectric function e e s, have been a useful tool for the interpretation of absorption 
spectral. Mean-field theories of the effective response function have also been proved to 
be useful to analyze electron energy-loss experiments in which swift electrons penetrate 
composites at random distances from the particles! - !. 

In the long- wavelength limit, the effective transverse dielectric function for a system of 
spherical particles was first derived by Maxwell-Garnett (MG)0, within a mean-field ap- 
proximation valid for small values of the volume occupied by the spheres. A generalization 
of the MG theory to finite wavelength, taking account of the spatial nonlocality in the 
response, has been recently presented by Barrera and Fuchsi. On the other hand, many 
attempts have been made to account, at large filling fractions, for higher multipole interac- 
tions, which are absent in the MG theory0-0. Recently, exact numerical calculations for 
periodic dielectric structures have been performed^, based in the long-wavelength limit of 
photonic band structure calculations, and the MG theory has been shown to offer a good 
approximation in the low-filling-ratio regime. Finally, the discovery of tubular fullerenes0, 
a few years ago, has opened a new focus of interest for the application of these theories to 
the case of cylindrical interfacesEiJii. 

In this paper we evaluate, on the basis of photonic band structure calculations, the 
long-wavelength limit of the effective transverse dielectric function of a composite made 
up of long metallic cylinders embedded in an otherwise homogeneous medium. We first 
demonstrate that the composite can be replaced by an effective homogeneous medium, and 
then we determine the transverse dielectric function of that medium from the electromagnetic 
modes supported by such a homogeneous system. We consider a spectral representation of 
the effective dielectric function, we present an evaluation of the parameters involved, and 
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we discuss the range of validity of a MG-like approach appropriate for cylindrical particles. 
We also address the impact of our results on absorption spectra and electron energy-loss 
phenomena. 



II. THEORY 

It was shown by Bergma ni and Miltonil that in the long-wavelength limit the effective 
dielectric function of any macroscopically homogeneous two-component system of dielectric 
constants e and eo and volume fractions / and 1 — / , respectively, can be expressed as a sum 
of simple poles which depend only on the microgeometry of the composite material and not 
on the dielectric functions of the components: 



e e fr — eo 



WE 



u — m. 



(1) 



where u is the spectral variable 



u = [1 - e / eo ] 



-i 



(2) 



m u are depolarization factors, and B v are the strengths of the corresponding normal modes, 
which all add up to unity: 



Similarly^, 



and 
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Ea = i. 
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In particular, if there is only one mode with a strength different from zero, one finds: 

1 



w- 



u — mi 



(6) 



and 

We consider a binary composite with a volume fraction 1 — / of insulator, with real 
dielectric constant eo, and / of a periodic system of long metallic cylinders of diameter d 
arranged in a square array with lattice constant a = xd (See Fig. 1). The diameter of 
the cylinders is taken to be small in comparison with the wavelength of the electromagnetic 
excitation and large enough that a macroscopic dielectric function e{u) is ascribable to the 
cylinders. For simplicity, the magnetic permeabilities are assumed to be equal to unity in 
both media. 

Now we take an electromagnetic (EM) wave normally incident on the structure, so that 
k x = k y = 0. For this propagation direction there are two different values of e e ff(u;) corre- 
sponding to different polarizations. In the case of EM waves polarized along the cylinders 
(s polarization), the presence of the interfaces does not modify the electric field, and one 
easily finds@ that the response of the composite is equivalent to that of a homogeneous 
medium having the effective transverse dielectric function of Eqs. (6) and (7) with mi = 
and rii = f. 

For EM waves polarized normal to the cylinders (p polarization), we first consider the 
case of a single cylinder embedded in an otherwise homogeneous medium. An elementary 
analysis of this classical problem shows that the electric field in the interior of the cylinder 

is! 

E in = ^^E out , (8) 

where E out represents the electric field intensity in the medium outside the cylinder. The 
effective dielectric function of isolated cylinders is then found to be given, again, by Eqs. 
(6) and (7), the depolarization factors now being m\ — n\ — 1/2. In the case of isolated 
particles only the dipole surface mode enters Eqs. (1) and (4) for this polarization. 

The interaction between the particles in a composite has been considered, over the years, 
within the well-known MG approximation. The basic assumption of this approach is that 



the average electric field within a particle located in a system of identical particles is related 
to the average field in the medium outside as in the case of a single isolated particle, thus 
only dipole interactions being taken into account. Hence, after replacing, for p polarization, 
E out of Eq. (8) by 

< E > -/E in 

^out - ' W 

< E > being the macroscopic effective electric field in the composite, the MG result for 
cylindrical inclusions is found, which can be represented by the use of Eqs. (6) and (7) with 

mi = i(l-/) (10) 



and 



11 = ^(1 + /)- (11) 



This is a special case of the generalized MG dielectric function derived in Ref. |25] for el- 
lipsoidal inclusions, which in the case of aligned spheroids distributed at random has been 
extended beyond the MG approximationil. Higher multipole interactions are neglected in 
these approaches. 

In order to compute, with full inclusion of higher multipoles, the effective dielectric func- 



tion of our periodic system, we first follow Ref. |2(J to calculate the photonic band structure. 
We fix the frequency, the local dielectric function inside the cylinders being therefore speci- 
fied, and we discretize the Maxwell equations within the unit cell. The discretized Maxwell 
equations provide a relationship between the electromagnetic fields on either side of the unit 
cell, i.e., z and z + a (see Fig. 1), and we exploit this to calculate, on applying Bloch's 
theorem, the eigenvalues of the so-called transfer matrix, which will give the band structure 
of the system. 

Homogeneous media are well known to support only two transverse modes. On the other 
hand, looking at the imaginary part of the eigenvalues of the system, we observe that over the 
entire range of frequencies there are, for each polarization, two degenerate EM Bloch waves 



with vectors k and —k for which k(u) roughly follows the dispersion relation of free light 
and which are dominant in the sense that they present the smallest attenuation (smallest 
imaginary part of k). Thus, we conclude that the composite can be replaced by an effective 
homogeneous medium, and we use the dispersion relation of these Bloch waves to compute 
the effective transverse dielectric function from 

k 2 c 2 

e cS (u) = —, (12) 

where c is the speed of light. 



III. RESULTS AND DISCUSSION 

We have taken the insulating component of the composite to have a real dielectric con- 
stant of eo = 1, and we have modeled the dielectric properties of the metal inside the 
cylinders by a simple free-electron (Drude) dielectric function: 

e("> ) = 1 " 7 - • s , (13) 

where u p is the bulk plasma frequency of the metal, and 7 represents an inverse electron re- 
laxation time. We have used the plasma frequency of the conduction electrons in Aluminum, 
hu p = 15.8eV, and 7 = 1.4eV. 

For s polarized electromagnetic waves, our numerical results for the effective dielectric 
function, as obtained from Eq. (12), accurately reproduce the exact results of Eq. (6) and 
(7) with m = and n = f, and this represents, therefore, a good check of our scheme. 

In the case of p polarized electromagnetic excitations, great care was exercised to ensure 
that our calculations are insensitive to the precise value of the number of mesh points in the 
unit celM The results presented below have been found to be well-converged for all filling 
fractions, and they all correspond to a sampling mesh of 180 x 1 x 1800. 

The optical absorption is directly given by the imaginary part of the effective dielectric 
function, Im [e e ff (a;)]. Thus, for p polarization the absorption spectra of isolated cylinders 
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(/ — > 0) exhibit a strong maximum, absent in the bulk metal, at the cylinder dipole resonance 
(e(u) + e = 0), i.e., at U\ = ^m\io v with m\ = 1/2. At higher concentrations of metal, 



Im[e e ff(u;)] still shows a dipole resonance at iO\ = Jm\{f)u p (mi(f) ^ 1/2), and a band of 



multipole resonances at u u = ^Jm v (f)u) p {y = 2,... ; m u > mi) is also formed broadened 
by electromagnetic interactions among individual cylinders. This is illustrated in Fig. 2, 
where our first principles calculations of Im[e e ff(u;)] // are exhibited, as a function of the 
frequency u, for various values of the ratio x between the lattice constant and the diameter 
of the cylinders. For x > 2 the effective dielectric function is well described by the MG 
approximation, and both our computed result and the MG approximation almost coincide 
for x = 6 with the isolated cylinder result, represented in Fig. 2 by a thin solid line. In Fig. 
3 our results are compared, for two different filling fractions, with the MG approximation, 
showing that together with the appearance of multipole resonances the depolarization factors 
mi(f) are, for x < 2, smaller than the MG depolarization factor of Eq. (10). 

In the case of a broad beam of charged particles penetrating the composite, the proba- 
bility per unit path length, per unit energy, for the projectile to transfer energy Two to the 
medium is^ 

2e 2 



h 2 v 2 



F{u), (14) 



F(u) = - dQ—^—lm \-e:} f (q, u)} , (15) 



where 



7r Jo Q 2 + ql 

v and e~jj-(q, to) being the velocity of the incident particles and the effective inverse longitu- 
dinal dielectric function of the composite, respectively. hQ represents the component of the 
momentum transfer, hq, located in a plane perpendicular to the incident beam direction, Q c 
is determined by q 2 c = Q 2 C + q%, q c representing the largest wave vector for which plasmons 
are well-defined excitations^, and 

q v = -■ 16 

v 
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In typical electron energy-loss experiments swift electrons of about lOOkeV penetrate 
the composite, and the main contribution to the integral of Eq. (15) comes, therefore, from 
very small values of the momentum transfer. Also, for small values of the adimensional 
parameter qa (qa < 1), a being the radius of the cylinder, the effective dielectric function 
of isolated cylinders has been proved to be rather insensitive to the introduction of a finite 
wave vectorial, and the energy-loss spectra of swift electrons is therefore expected to be, for 
small values of the radius of the cylinders (a < lnm)ill, well described by the q — > limit of 
the imaginary part of the effective inverse dielectric functionil. 

Eq. (7) with n = (dilute limit of n = /) and n = 1/2 exactly coincides, in the long- 
wavelength limit, with the result one obtains for the effective inverse longitudinal dielectric 
function of isolated cylinders when the momentum transfer is parallel to the axis and when 
it is located in a plane perpendicular to the axis of the cylinder, respectively^. In this limit 
(q — > 0) bulk plasmons are not excited, as a consequence of the so-called Bregenzung effect, 
and the pole in Eq. (7) at e(u) + e = 0, i.e., uj\ = ^JraxUp with m\ = 1/2, determines, for 
p polarization, the frequency of the so-called surface plasmon resonance at which electron 
energy-loss occurs. 

Fig. 4 shows, as a function of the frequency, our calculations of Im[— €^(u)]/f obtained 
for various filling fractions by varying the ratio x between the lattice constant and the diame- 
ter of the cylinders. In a homogeneous metal the so-called energy-loss function, Im[— e~ s (&)}, 
shows a single peak at the bulk plasmon resonance, and this peak persists in the compos- 
ite for x < 1, i.e., for filling fractions of metal for which the insulator no longer forms a 
connected medium. Besides this peak, there is also a band in Im[— e^(u})] from multipole 
contributions to the effective response at uj v = ^n v (f)uj p [y = 2, ... ; n„ < m). As in the 
case of Im [e e g(u;)], the energy- loss function is well described by the MG approximation in 
the low-filling-ratio regime, and our results begin to deviate from those obtained within the 
MG approximation at x ~ 2 (see Fig. 5). We also find that the depolarization factors n\(f) 
are larger than the MG depolarization factor of Eq. (11) for all volume fractions of metal 
for which the MG theory does not reproduce our results. 
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Our numerical calculations of the effective dielectric function, as obtained from Eq. (12), 
lead us to the conclusion that the depolarization factors m\ and ri\ of Eqs. (1) and (4), 
corresponding to dipolar modes, satisfy the following relation: 

ni = l-(D-l)mi, (17) 

where D represents the dimensionality of the inclusions, i.e., D = 3 in the case of spherical 
inclusions and D = 2 in the case of long circular cylinders. If all multipolar modes are 
neglected, then B\ = C\ = 1 in Eqs. (1) and (4) and a combination of these equations with 
Eq. (17) results, for D = 2, in the spectral representations of Eqs. (6) and (7) with the 
depolarization factors m\ and n\ given by Eqs. (10) and (11), i.e., the MG approximation. 
On the other hand, as long as multipolar modes give non-negligible contributions to the 
spectral representation of the effective response, i.e., B v ^ and C v ^ (u = 2,...), the 
strengths of the dipolar modes, B\ and Ci, become smaller than unity (see Eqs. (3) and 
(5)), and a combination of Eqs. (1) and (4) with Eq. (17) leads us to the conclusion that 
the depolarization factors rt\ and m 1 necessarily deviate from their MG counterparts. That 
a non-vanishing contribution from multipolar modes appears together with a deviation of 
the dipolar mode locations with respect to their MG counterparts is apparent from Figs. 3 
and 5. 

One of the main advantages of the spectral representation of Eqs. (1) and (4) is that the 
depolarization factors and the strengths of the corresponding normal modes depend only 
on the microgeometry of the composite material and not on the dielectric function of the 
components. On the other hand, the normal mode positions are determined by m v and n u , 
and we have thus represented in Fig. 6, in the case of p polarization, universal curves for 
the actual dipolar mode depolarization factors and strengths versus the filling fraction. It 
is obvious from this figure that the trend with increasing filling fraction is for the dipolar 
peaks in Im [e e e(a;)] / Im —€~£(u) to move from the cylinder dipole mode (mi — n\ — 1/2) 
to lower / higher energies. Notice that the MG results, also plotted in this figure, agree 
with our numerical results for / < 0.2 [x > 2). Nevertheless, at higher filling fractions the 
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actual depolarization factors approach more rapidly, as / is increased, the depolarization 
factors appropriate to the homogeneous metal (mi = and n\ — 1). In particular, for x < 1 
(/ > 7r/4) ni = 1, indicating that losses due to the bulk plasmon are expected when the 
metal forms a connected medium. 

Associated to the deviation of the actual dipolar mode positions from the MG results 
is the appearance of non-negligible multipolar mode strengths, B v [y = 2, ...) and C v [y = 
2,...), and, accordingly (see Eqs. (3) and (5)), a reduction in the dipolar mode strengths 
B\ and C\. These mode strengths determine the dipolar peak heights in Im [eefr(co')] and 



Im 



-e e ff K 



, represented in the inset of Fig. 6 as a function of /. For B\ = C\ 



all multipolar mode strengths being neglected, the dipolar peak heights in Im[e e ff(u;)] and 



Im 



Q 1 ' 



are found to be / H/ yjm{ and / H/ y/n^ } respectively, where H represents the 
peak height in the bulk energy-loss function (in the case of the Drude dielectric function 
of Eq. (13), H = ujp/j). The actual dipolar peak heights divided by f H are represented 
by stars and dots, together with the B x — C\ — \ approximation, i.e., 1/ \fm{ and \j^fn{, 
showing, therefore, the magnitude of the actual dipolar mode strengths involved in the 
spectral representations of Eqs. (1) and (4). 

In summary, we have evaluated, on the basis of photonic band structure calculations, the 
effective long-wavelength dielectric response function of a system of long metallic cylinders. 
We have considered a spectral representation of both direct and inverse effective dielectric 
functions, of interest in the interpretation of absorption spectra and electron energy-loss 
experiments, respectively. We have analyzed the surface mode strengths and positions, and 
the effect of the electromagnetic interactions between the cylinders has been investigated. 
We have concluded that MG results are good as long as the distance between the axis of 
neighbouring cylinders is larger than twice the diameter of the cylinders. This is in contrast 
with the results we have also obtained for long dielectric cylinders. In this case we have 
found the MG results to be good for almost all the filling ratios for which the cylinders are 
not touching, as happens in the case of a periodic array of dielectric spheres@. 
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FIGURES 

FIG. 1. Periodic array of metallic cylinders of diameter d arranged in a square array with lattice 
constant a. The cylinders are infinitely long in the y-direction. 

FIG. 2. The imaginary part of the effective long-wavelength dielectric function of the periodic 
system described in Fig. 1, for p polarized electromagnetic excitations. Full thin line: the isolated 
cylinder result (/ — >• 0). Dotted, long-dashed, short-dashed, dashed-dotted and thick full lines: our 
numerical results for volume filling fractions of 2.2%, 19.6%, 54.5%, 68.6% and 74.0%, respectively. 

FIG. 3. The same as in Fig. 2, for filling fractions of 54.5% and 74.0%. Full thin and thick 
lines: our numerical results. Dotted and dashed lines: MG results. 

FIG. 4. The effective energy-loss function of the periodic system described in Fig. 1, for p 
polarized electromagnetic waves. Full thin line: the isolated cylinder result (/ — ► 0). Dotted, 
long-dashed, short-dashed, dashed-dotted and thick full lines: our numerical results for volume fill- 
ing fractions of 2.2%, 19.6%, 54.5%, 68.6% and 78.5%, respectively. The bulk energy-loss function 
is represented bay a dotted line. 

FIG. 5. The same as in Fig. 4, for filling fractions of 54.5% and 78.5%. Full thin and thick 
lines: our numerical results. Dotted and dashed lines: MG results. 

FIG. 6. Dipolar mode positions (depolarization factors), as a function of the volume filling 
fraction /. Solid lines: our numerical results for mi (the line below the dilute limit of 1/2) and 
n\ (the line above the dilute limit of 1/2). Dotted lines represent the MG depolarization factors 
given by Eqs. (10) and (11). Dipolar mode strengths are represented in the inset. Dots and stars: 
our numerical results for the dipolar peak heights in Im[e e ff(u;)] /(/ H) and Im 
respectively. Full thick and thin lines: l/y/m 1 and l/y / n 1 , respectively. Here, H = Up/j. 
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